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ABSTRACT 

We present an analysis of the X-ray bright point source population in 43 massive clusters 
of galaxies observed with the Chandra X-ray Observatory. We have constructed a catalog 
of 4210 rigorously selected X-ray point sources in these fields, which span a survey area 
of 4.2 deg . This catalog reveals a clear excess of sources when compared to deep blank- 
field surveys, which amounts to roughly 1 ± 0.4 additional source per cluster, likely Active 
Galactic Nuclei (AGN) associated with the clusters. The excess sources are located within the 
overdensity radius r^oo, with the largest excess observed near the cluster centers. The average 
radial profile of the excess X-ray sources of the cluster are well described by a power law 

(N(r) ~ r@) with an index of /3 0.6. An initial analysis using literature results on the mean 

profile of member galaxies in massive X-ray selected clusters indicates that the fraction of 
galaxies hosting X-ray AGN rises with increasing clustercentric radius, being approximately 
3 times higher near r 5 oo than in the central regions. This trend is qualitatively similar to that 
observed for star formation in cluster member galaxies. 



1 INTRODUCTION 



A critical element of galaxy evolution models is the influence of the 
local environment. Environmental effects are expected to be partic- 
ularly pronounced in galaxy clusters, where hundreds or thousands 
of galaxies are in close proximity to one another, and are embedded 
in a hot, diffuse intracluster medium (ICM). Previous studies have 
shown that the ICM has clear effects on member galaxy properties: 
the star formation rate of galaxies in clusters is lower when com- 
pared to field galaxies, resulting in cluster member s being observ- 
ably redder and more elliptical (e.g. |PressleJll980h . This suggests 
that the ICM is efficient at stripping cold gas from galaxies and 
suppr essing star formation, at least since z ~ 1 (e.g. lElbaz et all 
120071) , 

Theoretical and observational studies have shown that the evo- 
lution of Active Galactic N uclei (AGN) and th eir host galaxies are 
linked to one another (e.g. ISilk & Reeslll998h . The masses of the 
supermassive black holes at the centers of nearby galaxies correlate 
tightl y with the masse s of the i r galaxy bulges (e.g.lMagorrian et al.l 
1 19981 : iGebhardt et all l2000l : iMarconi & Hunl 120031) suggesting 
a connection between their growth histories. Feedback between 
AGN and their host galaxies has also been observed to play 
a key role in suppressing star formation (e.g. Kauffmann et al.l 



dio studies have led to th e conclusion tha t there are at least two 
distinct AGN modes (e . g. iBest et alj|2007 



007) : 

iMerloni & Heinj|2008l : TBest & Heckmanll2012l : lFabianll2012l) : a 
radiatively efficien t mode associated wi t h X-ray and optically- 
selected AGN (e.g lAlexander et al.ll2003l: iKauffmann e7aLll2003l: 
iBauer et al. 1 1200 % IXue et al.1 1201 01) and a radiatively inefficient 
mode, associated with relatively strong jets and radio emission 



Hardcastle et~ai1l2007l: 



halos (e.g. Allen et alj|200d: 


Balmaverde et alj|2008l: iDunn et al.l 


l20ld:lBest & Heckmanl 20121: 


Pellegrini et alj|2012l). While the ra- 



a Key role m su ppressing star tormatton (e.g. Kauttmann et al. 
200a:lHopldnd2004lCro"ton et alj200fj;lHopkins & Beacoml2006l : 



McNamara & Nulsenl 120071 : lFabianll2012l) . X-ray, optical, and ra- 



diatively efficient mode appears physically related to the accretion 
of cold gas onto the black hole, the radiatively inefficient mode may 
be associated with the accretion of hot X-ray emitting gas. The ki- 
netic energy injected into the surrounding environment by the jets 
emitted in this radiatively inefficient mode is, in principle, suffi- 
cient to prevent the cooling of the hot gas and suppress subsequent 
star formation (e. g.lAilen et alj|2006l : iMcNamara & Nulsenll2007l: 
IDunn et aljirJTol : |Fabianll20 1 2h . 

The differences between the occurrence of AGN in clus- 
ter galaxies versus the field have not been well established . 
Early studies of cluster membe r galaxies (e.g. IOsterbrock|[l960] : 
lGislerll978l : lDressler et al.ll985f) determined that optically selected 
AGN were relatively rare i n clusters. More recent work utiliz- 
ing large area surveys (e.g. IBest et ai] |200l l Arnol d et al. I l2009l: 



© 0000 RAS 



iKoulouridis & Plionisll201ol : I von der Linden et"al1l2010l) have ex- 
panded these findings, and prompted vigorous debate. Some stud- 
ies have argued that optically luminous AGN (Lroim > 10 7 ^n ) 
reside predominantly outside of clusters jKauffmann et alj|2004h . 
while others have suggested th at the optical AGN activity is in- 
dependent of the environment dMiller et alj l2003h . Observations 
at radio wavelengths suggest that cluster environments are more 
conducive to AGN activity. Brightest Cluster Galaxies (BCG's), in 
particular, appear to be likely to host radio-luminous AGN while 
less likely to host an optically lumin ous AGN when compared to 
field galaxies of the same mass (e.g.lBestll20q4lBest et al.ll2007l; 
von der Linden et al.1 120071 : iDunn et alj |2010| : I Best & Heckmar] 



20121) . 



Comparisons between the cluster and field populations of 
AGN using X-ray observations have led to differing conclusions, 
and have primarily investigated whether or not X-ray point sources 
are more frequently o bserved in the vicinity of galaxy clusters. 
ICappelluti et alj J2005h identify an excess of point sources in four 
separate cluster fields relative to expectations from field surveys, 
although for six other cluster fields in their study no such excess 
was seen. Excesses have also been observed i n particular cluste r 
fields such as 3C29 5 RX J003033 2+26 1819 l lCappi et alj|200lh . 
and MS 1054-0321 Johnson et ai1l2003h . but no excess of X-ray 
point sources was detected in the vicinit y of galaxy clus ters ob- 
served in the large field ChaMP survey dKim et alj|2004h . Other 
studies using larger samples of galaxy clusters have argued for a 
statistical excess in their fields, and have shown that the excess of 
sources is primarily located within the central 1-2 Mpc o f each 
cluster (e.g. lRuderman & Ebeling|2005l ; lGilmour et alj20"oH) . Mul- 
tiwavelength studies of cluster galaxies incorporating optical spec- 
troscopy along with X-ray analysis have suggested that the fraction 
of galaxies hosting X-ray bright AGN in cluste rs may increase sub- 
stanti ally with redshift, by a factor of ~ 8 (e.g. lMartini et alj|2007l . 
120091). an effect that was also observed in the large Bootes field 
dGalametz etai]|2009l) . 

X-ray surveys of AGN offer hi gher source densities of 
AGN than other wavelengths (see e.g. IBrandt & Hasingerl 120051 ; 



iBrandt & Alexander] 1 20 lof for reviews). Deep field studies of the 
X-ray sky such as the Chandra Deep Field North and Chan- 
dra De ep Field South (he r eafter CDFN and CDFS, respectively 
) (e.g IBrandt et all I2OO1I; iGiacconi et al.1 120021 : lAlexander et"al1 
I2OO3I : iLuo et al.ll2008l ; IXue et alj 2011) have precisely measured 



the surface density of X-ray point sources down to fluxes of ~ 
3 x 10~ 17 ergcirT 2 s~' in the 0.5 - 8.0 keV band. Above fluxes of 
~ 10~ 14 ergcirT 2 s -1 in the 0.5 - 8.0keV band, the majority of 
these X-ray point sources are unobscured or mildly obscured AGN 
(e.g. iBade et al.|[l998l ; ISchmidt et afll 19981 ; lAkivama et aT]|20031) . 
At lower fluxes, however, other source classes such as obscured 
AGN, sta rburst galaxies, and normal galaxies bec ome more promi- 
nent (e.g. iBauer et al]|2004 iLehmer et all 1201 2h . Although X-ray 
surveys have measured the largest densities of reliably identified 
AGN at any wavelength, multiwavelength follow-up is necessary 
to test fully the impact of the cluster environment on the evolution 
of AGN. Field sources cannot be easily distinguished from sources 
intrinsic to the cluster using X-ray observations alone, even when 
incorporating X-ray spectroscopy. The typical spectrum for unob- 
scured AGN in the ~ 2 - lOkeV energy band is characterized, 
at least to first order, by a power-law continuum, N(E) ~ E~ r , 
where the photon index r is typically in the rang e of ~ 1.5 - 2.5 
jNandra & Pounds! 1 1994 iReeves & Turner! l200d ; IShemmer etall 
l2005l : |just et al.ll2007t) . making it difficult to identify spectral fea- 
tures and so measure source redshifts. X-ray observations also pro- 



vide little insight into the underlying population of galaxies in a 
given cluster, which diminishes the physical significance of ob- 
serving a larger number of AGN in these fields. A multiwavelength 
study of galaxy clusters utilizing both X-ray and optical observa- 
tions is well suited to offer a more complete picture of the AGN 
and galaxy populations within the cluster. 

Here, we present the first results of a study designed to select 
and observe AGN and galaxies in clusters at both X-ray and op- 
tical wavelengths. Particularly, we utilize deep observations with 
the Chandra X-ray Observatory and the Subaru optical telescope 
of forty-three of the most massive galaxy clusters in the redshift 
range 0.2 < z < 0.7, which also have deep, multi-filter Subaru 



optical imaging (see 


von der Linden et alj|2012l 


Applegate et al. 2012 


). An important element of 



of the clusters have well-determined masses based on X-ray mass 
proxies, allowing us to examine the AGN and galaxy properties as 
a function of overdensity in the cluster (i.e. relative to the virial ra- 
dius as opposed to a metric aperture). In this first paper, we focus 
on the spatial distribution of X-ray point sources within the cluster 
fields. 

The structure of this paper is as follows: Section 2 discusses 
the sample of clusters utilized in this study. Section 3 discusses the 
initial processing of the X-ray data, while Section 4 discusses the 
production of X-ray point source catalogs for each cluster field. 
Section 5 addresses the calculation of the point-source sensitivity 
map for each cluster field. Section 6 presents our results on the cu- 
mulative number counts and radial profiles of X-ray point sources 
in these fields. In Section 7 we discuss the implications of these 
results on the evolution of X-ray selected AGN in galaxy clusters. 
For calculating distances, we assume a ACDM cosmological model 
with n m =0.3, n A =0.7, and H = 70kms-' Mpc" 1 . 



2 THE CLUSTER SAMPLE 

The clusters in our study have been drawn from three wide-area 
cluster surveys d erived from the ROSAT All Sky Survey (RASS; 
TruempeJ 1 19931): the ROSAT Brightest Cluster Sample (BCS; 
Ebeling etalJI 19981); the ROSAT-E SO Flux-Limited X-ray Sample 
(REFLEX: lB6hringeret"ai1l2004h: and the MAssive Cluster Sur- 
vey (MACS: lEbeling et alj|2007l.l20ld) . Each sample covers a dis- 
tinct volume of the Universe: the BCS covers the northern sky at 
z < 0.3; REFLEX covers the southern sky at z < 0.3; and MACS 
covers higher redshifts, 0.3 < z < 0.7, at declinations > 40°. Clus- 
ters included in these samples have been instrumental in the most 
recent cosm ological studies u tilizing growth of structure tests (for 
a review see lAllen etalj|201ll) . 

For this paper, we restrict ourselves to the subset of these 
clusters with the most thorough multiwavel ength follow-up me a- 
surements, including Chandra observations dMantz et alj|2010dlbh 
and deep multi-filter optical imaging with the Subaru telescope. 
In detail, we require at least 10 ks of exposure time with Chan- 
dra; robust measurements of the cluster mass and virial radius 
from iMantz et~aL I feOlOdlbh ; and that the clusters are included 
in t he weak-lensing analysis of von d er Linde n et al. I d2012l) (see 
also lKellv et alJl2012l ; lApplegate et alj|2012l) l'l. In total. 43 unique 
galaxy clusters have been analyzed here, with redshifts ranging 



1 As discussed in these papers, the clusters chosen for optical follow-up 
were selected based on their X-ray properties. This sample is therefore ex- 
pected to be unbiased with respect to the optical properties of the clusters. 



from 0.2 < z < 0.7. General information for the clusters and the 
Chandra data sets used may be found in Table[T] 

The clusters studied here are among the most massive and X- 
ray luminous clusters in the Universe, and host high densities of 
both galaxies and ICM. We therefore expect the influences of the 
cluster environment to be pronounced in this sample. With mea- 
surements of 7-500 for each cluster we are able to relate observed 
trends in the AGN population to the virial radii of the clusters. Mass 
measurements and th e associated overden sity radii, r 500 , for each 
cluster are taken from lMantz etafl d2010dlbh .Fl The typical uncer- 
tainties in measurements of r 50 o are of order ~ 10%. The r 500 values 
and X-ray centroids for the clusters are summarized in Table[T] 



3 Chandra DATA REDUCTION 

All of the galaxy cluster fields were observed with the Advanced 
CCD Imaging Spectrometer (ACIS) aboard Chandra. The major- 
ity of the observations utilize the ACIS-I chip array, but we also 
include observations utilizing ACIS-S. For ACIS-S observations, 
we only utilize chips SI, S2, and S3, for which the most accurate 
calibration products are available. The standard level- 1 event lists 
produced by the Chandra pipeline processing were reprocessed us- 
ing the CIAO (version 4.3) software package, including the appro- 
priate gain maps and calibration products (CALDB version 4.4.6). 
Bad pixels were removed and standard grade selections were ap- 
plied. The data were cleaned to remove periods of anomalously 
high background, using the standard energy ranges and binning 
methods recommended by the Chandra X-ray Center. The net ex- 
posure times after processing are summarized in Table Q] Images 
for each cluster observation were produced in each of three energy 
bands, hereafter designated the full band (0.5 - 8.0 keV), the soft 
band (0.5 - 2.0 keV), and the hard band (2.0 - 8.0 keV). 

For each observation, an exposure map was produced in each 
energy band, calculating the product of quantum efficiency and ef- 
fective area across each field of view. Each exposure map was spec- 
trally weighted over its respective energy band, assuming a pho- 
ton index of T = 1.4 and a hydrogen column density fixed to 
the value measured in the Leiden-Argentine-Bonn survey (LAB, 
iKalberla et al .120051) . Effective exposure time maps were calculated 
for each observation and energy band by dividing the exposure map 
by its maximum effective area and multipl ying by the exposure 
time, similar to the procedure discussed in iHornschemeier et all 
d200ll) . In instances where there was more than one Chandra ob- 
servation of a cluster, images and exposure maps were created for 
each observation individually, and then reprojected into a common 
aspect solution and combined into a single coadded cluster im- 
age/exposure map. These coadded images and exposure maps only 
include active CCD chips from the primary field of view for each 
observation (i.e. chips 10-13 for ACIS-I observations, and chips S 1- 
S3 for ACIS-S observations). 

2 The overdensity radius r'A is denned as the radius where the enclosed 
average mass density is equal to A times the critical density of the universe 
at the cluster's respective redshift, p c (z). We define the virial radius for this 
study as rioo- The corresponding mass is defined as Ma = 4/3nAp c (z) 
ta 3 . The mass range extends from 4 X 10 14 M© < M500 < 2 X 10 I5 M G and 
the overdensity radii range from 1.06 Mpc < /'500 < 1.74 Mpc. 

3 For this study, the canonical AGN source was chosen to have a photon 
index of V = 1.4, to remain consistent with other studies of AGN popula- 
tions in the field, in particular the CDFS (e.g. lXue et alj201 U lLehmer et alj 
l20Tj . 



4 CONSTRUCTING THE X-RAY POINT SOURCE 
CATALOG 

4.1 Initial Candidate X-Ray Point Sources 

An initial catalog of candidate source s was produced for ea ch clus- 
ter using the WAVDETECT algorithm dFreeman et"ai1l2002l) and the 
full-band (0.5 - 8.0 keV) image. Our WAVDETECT runs utilized 
a "V2 sequence" of wavelet scales (i.e. 1, V2, 2, 2V2, 4, 4V2, 
8, 8 V2, and 16 pixels) and a false-positive probability threshold 
of 10~ 5 . This false-positive threshold is more liberal in including 
faint sources than that commonly used in other studies and allows 
for recovering legitimate sources that fall below form ally stricter 
thresholds (e.g. lAlexander et al.ll200ll ; rXue et alfcoill) . The trade- 
off, however, is that we incur an appreciable number of spurious 
sources in our initial catalogs, many of which are associated with 
structures in the cluster ICM. Where all observations of a given 
cluster utilized a common detector (e.g. all observations were per- 
formed with ACIS-I), we utilized WAVDETECT's PSF information 
in determining the positions and size of candidate sources. The 
WAVDETECT software employs a single PSF model, and thus reli- 
able PSF information cannot be supplied into WAVDETECT for the 
seven clusters with mixed ACIS-I and ACIS-S observations. Can- 
didate AGN for these clusters were therefore identified assuming 
the smallest possible wavelet scale (1 pixel), resulting in a boosted 
number of candidate sources. However, our subsequent screening 
processes are efficient at removing spurious sources, leading to the 
results for these clusters being consistent with the rest. Our initial 
WAVDETECT catalogs include a total of 6266 sources across all 
cluster fields. Tests show that it is unnecessary to subsequently run 
WAVDETECT on the hard-band and soft-band images separately, 
as only a few additional sources per cluster field are identified in 
these images. Moreover, the majority of these additional sources 
detected in the hard or soft bands are subsequently rejected after 
further analysis. 

4.2 Improving the Initial X-ray Point Source Catalogs 

Given our liberal false-positive probability threshold, incomplete 
incorporation of the Chandra PSF, and the bright background ICM 
present at the center of every field of view, we expect that the ini- 
tial catalogs produced by WAVDETECT have an appreciable number 
of spurious sources. Further analysis is required to ensure source 
validity. To achieve this, every candidate source in the initial cat- 
alog is re-analyzed on an observation-by-obse rvation basis usin g 
the a? point-source analysis software package^ dBroos et"aT]|2010h . 
This offers several key improvements, specifically: (1) utilization 
of source and background regions that approximate the shape of 
the PSF calculated from ray-tracing simulations, accommodating 
neighboring sources and CCD chip boundaries as needed; (2) a 
multi-stage approach to source detection that refines the source cat- 
alog based on each source's binomial no-source probability (i.e. the 
probability of a source not being real given local backgrounds); and 
(3) better source-position determinations that improve the accuracy 
of each source's as trometry and photometry. We refer the reader to 
iBroos et al.1 d20ld) for documentation about a?, although the most 
relevant aspects of the analysis are discussed below. 

Our analysis pipeline has three main stages and follows the 

4 The ACIS Extract software package and User's Guide are available at 
://www.astro.psu.edu/xray/acis/acis_analysis.html 



http 



Table 1. Summary of the cluster sample and Chandra data. A total of 43 clusters of galaxies and 89 independent Chandra observations were used. The 
columns are (1) cluster name; (2) redshift; (3) & (4) right ascension and declination of the cluster X-ray centroid, in degrees (J2000); (5) the Chandra 
observation number (OBS ID#); (6) observation date; (7) primary detector array used; (8) exposure time in ks, after cleaning; (9) the measured value of the 
cluster mass M500, in units of 10 14 M o ; (10) measured value of r$QQ for the cluster, measured in Mpc; and (11) the equivalent Galactic hydrogen column 
density in the direction of the cluster, in units of 10 20 atoms cirT 2 . 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Cluster Name 


Z 


RA 


DEC 


OBS ID # 


Obs Date 


Detector 


Exposure (ks) 


M 50 o(10 14 M G ) 


r5oo(Mpc) 


Nn (10 20 cirT 2 ) 


Abell 209 


0.206 


22.971 


-13.613 


3579 


2003-08-03 


ACIS-I 


9.98 


12.62 


1.531 


1.44 










522 


2000-09-09 


ACIS-I 


9.96 








Abell 963 


0.206 


154.264 


39.047 


7704 


2007-02-18 


ACIS-I 


5.06 


6.82 


1.247 


1.25 










903 


2000-10-11 


ACIS-S 


36.28 








Abell 2261 


0.224 


260.612 


32.132 


5007 


2004-01-14 


ACIS-I 


24.31 


14.42 


1.588 


3.19 










550 


1999-12-11 


ACIS-I 


9.06 








Abell 2219 


0.228 


250.084 


46.708 


896 


2000-03-31 


ACIS-S 


42.29 


18.87 


1.737 


1.76 


Abell 2390 


0.233 


328.404 


17.695 


4193 


2003-09-11 


ACIS-S 


95.06 


15.15 


1.613 


6.21 










500 


2000-10-08 


ACIS-S 


9.82 
















501 


1999-11-05 


ACIS-S 


9.04 








RXJ2129.6+0005 


0.235 


322.415 


0.088 


552 


2000-10-21 


ACIS-I 


9.96 


7.71 


1.285 


3.63 










9370 


2009-04-03 


ACIS-I 


29.64 








Abell 521 


0.248 


73.530 


-10.223 


430 


2000-10-13 


ACIS-S 


39.11 


11.43 


1.459 


4.78 










901 


1999-12-23 


ACIS-I 


38.63 








Abell 1835 


0.253 


210.258 


2.877 


6880 


2006-08-25 


ACIS-I 


117.91 


12.29 


1.493 


2.04 










6881 


2005-12-07 


ACIS-I 


36.28 
















7370 


2006-07-24 


ACIS-I 


39.50 








Abell 68 


0.255 


9.274 


9.160 


3250 


2002-09-07 


ACIS-I 


9.98 


7.61 


1.270 


4.96 


Abell 697 


0.282 


130.740 


36.365 


4217 


2002-12-15 


ACIS-I 


19.51 


17.10 


1.647 


2.93 


Abell 2537 


0.297 


347.091 


-2.191 


4962 


2004-09-09 


ACIS-S 


36.19 


7.15 


1.225 


4.62 










9372 


2008-08-11 


ACIS-I 


38.50 








MS2137.3-2353 


0.313 


325.063 


-23.661 


4974 


2003-11-13 


ACIS-S 


57.38 


4.74 


1.063 


3.76 










5250 


2003-11-18 


ACIS-S 


40.54 
















928 


1999-11-18 


ACIS-S 


43.60 








MACSJ1931.8-2634 


0.352 


292.956 


-26.576 


3282 


2002-10-20 


ACIS-I 


13.59 


9.94 


1.339 


8.31 










9382 


2008-08-21 


ACIS-I 


98.92 








MACSJ1 115.8+0129 


0.355 


168.966 


1.498 


3275 


2003-01-23 


ACIS-I 


15.90 


8.64 


1.278 


4.34 










9375 


2008-02-03 


ACIS-I 


39.62 








RXJ1532.9+3021 


0.363 


233.224 


30.349 


1649 


2001-08-26 


ACIS-S 


9.36 


9.48 


1.312 


2.30 










1665 


2001-09-06 


ACIS-I 


9.97 








ZWCL1953 


0.378 


132.529 


36.072 


1659 


2000-10-22 


ACIS-I 


24.86 


10.16 


1.336 


2.96 










7716 


2006-12-20 


ACIS-I 


6.98 








MACSJ0949.8+1708 


0.384 


147.465 


17.118 


3274 


2002-11-06 


ACIS-I 


14.31 


11.35 


1.380 


3.08 


MACSJ1720.2+3536 


0.387 


260.069 


35.606 


3280 


2002-11-03 


ACIS-I 


20.84 


6.31 


1.135 


3.46 










6107 


2005-1 1-22 


ACIS-I 


33.88 








MACSJ 173 1.6+2252 


0.389 


262.913 


22.863 


3281 


2002-11-03 


ACIS-I 


20.50 


12.82 


1.427 


4.99 


MACS J22 11.7-0349 


0.396 


332.941 


-3.828 


3284 


2002-10-08 


ACIS-I 


17.73 


18.06 


1.608 


5.53 


MACSJ0429.6-0253 


0.399 


67.400 


-2.884 


3271 


2002-02-07 


ACIS-I 


23.16 


5.76 


1.097 


4.34 


MACSJ2228.5+2036 


0.411 


337.136 


20.620 


3285 


2003-01-22 


ACIS-I 


19.85 


14.67 


1.491 


4.26 


MACSJ045 1.9+0006 


0.429 


72.977 


0.105 


5815 


2005-01-08 


ACIS-I 


10.21 


6.33 


1.118 


6.85 


MACSJ1206.2-0847 


0.439 


181.551 


-8.801 


3277 


2002-12-15 


ACIS-I 


23.45 


19.16 


1.612 


4.35 


MACSJ0417.5-1154 


0.443 


64.393 


-11.907 


11759 


2009-10-28 


ACIS-I 


51.35 


22.06 


1.689 


3.31 










12010 


2009-10-29 


ACIS-I 


25.78 
















3270 


2002-03-10 


ACIS-I 


12.01 








MACSJ2243.3-0935 


0.447 


340.839 


-9.595 


3260 


2002-12-23 


ACIS-I 


20.50 


17.35 


1.555 


4.02 



recent study of the 4-Ms CDFS dXue et al.ll201lh . modified to ac- 
commodate the higher background rates and shorter exposure times 
of our cluster observations. Unless otherwise noted, all calculations 
were performed in the full band (0.5 - 8.0 keV). 



4.2.1 Removing False Sources 

For each observation of every candidate point source, ae is used to 
compute a polygonal source-extraction region based on models of 
the Chandra PSF at the source position. These models use sim- 
ulations of the Chandra High Resolution Mirror Assembly using 
the MARX ray-tracing simulator. Polygonal extraction regions were 
constructed that approximate the ~ 90% encircled energy fraction 
(EEF) contour of the local PSF at 1.497 keV. When dealing with 



Table 1. Continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


Cluster Name 




RA 


DEC 


OBS ID # 


Obs Date 


Detector 


Exposure (ks) 


M 50 o(10 14 M s ) 


r 500 (Mpc) 


Nh (10 20 cm~ 2 ) 


MACSJ0329.6-0211 


0.450 


52.422 


-2.195 


3257 


2001-11-25 


ACIS-I 


9.86 


7.89 


1.194 


4.64 










3582 


2002-12-24 


ACIS-I 


19.84 
















6108 


2004-12-06 


ACIS-I 


39.64 
















7719 


2006-12-03 


ACIS-I 


7.08 








RXJ1347.5-1145 


0.451 


206.878 


-11.752 


3592 


2003-09-03 


ACIS-I 


57.71 


21.71 


1.674 


4.60 










506 


2000-03-05 


ACIS-S 


8.93 
















507 


2000-04-29 


ACIS-S 


9.99 








MACSJ1621.3+3810 


0.463 


245.353 


38.169 


10785 


2008-10-18 


ACIS-I 


29.75 


5.89 


1.078 


1.13 










3254 


2002-10-18 


ACIS-I 


9.84 








... 








6109 


2004-12-11 


ACIS-I 


37.54 


... 














6172 


2004-12-25 


ACIS-I 


29.75 
















9379 


2008-10-17 


ACIS-I 


29.91 








MACSJ1 108.8+0906 


0.466 


167.229 


9.100 


3252 


2002-11-17 


ACIS-I 


9.94 


7.73 


1.179 


2.22 










5009 


2004-02-20 


ACIS-I 


24.46 








MACSJ1427.2+4407 


0.487 


216.816 


44.125 


6112 


2005-02-12 


ACIS-I 


9.38 


6.35 


1.095 


1.19 










9380 


2008-01-14 


ACIS-I 


25.81 
















9808 


2008-01-15 


ACIS-I 


14.93 








MACSJ2214.9-1359 


0.502 


333.738 


-14.003 


3259 


2002-12-22 


ACIS-I 


19.47 


13.16 


1.387 


2.88 










5011 


2003-11-17 


ACIS-I 


18.52 








MACSJ091 1.2+1746 


0.505 


137.795 


17.775 


3587 


2003-02-23 


ACIS-I 


17.87 


8.96 


1.220 


3.28 










5012 


2004-03-08 


ACIS-I 


23.79 








MACSJ0257. 1-2325 


0.505 


44.287 


-23.434 


1654 


2000-10-03 


ACIS-I 


19.84 


8.51 


1.198 


2.08 










3581 


2003-08-23 


ACIS-I 


18.47 








MACSJ0454. 1-0300 


0.538 


73.547 


-3.014 


529 


2000-01-14 


ACIS-I 


13.90 


11.46 


1.308 


3.92 










902 


2000-10-08 


ACIS-S 


44.19 








MACSJ1423.8+2404 


0.543 


215.949 


24.078 


1657 


2001-06-01 


ACIS-I 


18.52 


6.65 


1.088 


2.20 










4195 


2003-08-18 


ACIS-S 


115.57 








MACSJ1 149.5+2223 


0.544 


177.397 


22.401 


1656 


2001-06-01 


ACIS-I 


18.52 


18.66 


1.533 


1.92 










3589 


2003-02-07 


ACIS-I 


20.04 








MACSJ0717.5+3745 


0.546 


109.383 


37.755 


1655 


2001-01-29 


ACIS-I 


19.87 


24.87 


1.688 


6.64 










4200 


2003-01-08 


ACIS-I 


59.16 








MS0015.9+1609 


0.547 


4.639 


16.436 


520 


2000-08-18 


ACIS-I 


67.41 


16.47 


1.469 


3.99 


MACSJ0025.4-1222 


0.585 


6.374 


-12.379 


10413 


2008-10-16 


ACIS-I 


75.63 


7.59 


1.119 


2.50 










10786 


2008-10-18 


ACIS-I 


14.12 
















10797 


2008-10-21 


ACIS-I 


23.85 
















3251 


2002-11-11 


ACIS-I 


19.32 
















5010 


2004-08-09 


ACIS-I 


24.82 








MACSJ2129.4-0741 


0.588 


322.357 


-7.691 


3199 


2002-12-23 


ACIS-I 


19.85 


10.65 


1.250 


4.33 










3595 


2003-10-18 


ACIS-I 


19.87 








MACSJ0647.7+7015 


0.592 


101.957 


70.248 


3196 


2002-10-31 


ACIS-I 


19.27 


10.88 


1.257 


5.40 










3584 


2003-10-07 


ACIS-I 


19.99 








MACSJ0744.8+3927 


0.698 


116.217 


39.457 


3197 


2001-11-12 


ACIS-I 


20.23 


12.53 


1.264 


5.66 










3585 


2003-01-04 


ACIS-I 


19.85 
















6111 


2004-12-03 


ACIS-I 


49.50 









crowded sources that have overlapping regions, a; utilizes smaller 
extraction regions (corresponding to ~ 40 - 75% EEFs) chosen 
to be as large as possible without overlapping. These source re- 
gions were used throughout the analysis, including the determina- 
tion of the source's local background. For our analysis, the as "BET- 
TER_B AC KGROUNDS ' ' algorithm was utilized for background 
extraction, which computes the background counts within external 
background regions surrounding the source, accounting for con- 
tributions from the source and neighboring sources. These back- 
ground regions were required to contain a minimum of 100 counts, 
and had a median of 1 1 1 counts. In order to provide the most accu- 
rate possible backgrounds, front-illuminated and back-illuminated 



chips in a particular observation were treated as separate observa- 
tions. 

In order to refine the catalog and remove spurious sources, we 
utilize the no-source binomial probability to determine the likeli- 
hood that source c ounts are due to fluctua tions in the background 
(see Appendix A o f lWeisskopf et alj|2007t) . If there are a total of S 
counts in a source region (subtending a solid angle of Q slc ) and 6 cxt 
counts in the external background region (subtending a solid angle 
of £2 cxt ), the binomial probability can be calculated as: 



I 



AM 



(N-X)lXl 



P x d-Pf- X 



(1) 



In this equation, N is the total number of counts in the 



C 
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Figure 1. The distribution of source detection probabilities and fluxes, (a): 
Distribution of X-ray no-source binomial probabilities, Pj,, for all sources 
before inclusion in the final catalog. For this figure, probabilities were de- 
termined from counts in the full band (0.5 - 8.0 keV). The minimum thresh- 
old chosen for inclusion in the final catalog, Pb < 10~ 3 , is denoted by the 
dashed vertical line. For presentation purposes, we set log Pb = -20 for all 
sources where log Pb ^ -20. 



source+background region (S+B exl ), and p= 1/(1 + BACKSCAL) is 
the probability that a background count is located within the source 
region (thus contributing to S), where BACKSCAL=C1 M /Q. src . 
Since a? calculates the expected number of background counts in 
the source region (6 src ) from this same area ratio, BACKSCAL 
is also equivalent to the ratio of background counts in the exter- 
nal background region and source regions (B exl /B slc ). The median 
value of BACKSCAL for all sources included in the final catalog is 
roughly 62. In the first stage of the pipeline, this probability was 
determined for each source and all sources with P b > 0.01 (i.e. a 
probability of being a background fluctuation of more than one per- 
cent) were removed. This ensures that the source and background 
regions for "real" sources are not influenced by the presence of 
spurious sources nearby. Although the results vary from cluster to 
cluster, a total of 1407 sources (22%) of the initial WAVDETECT 
sources were removed by this step. Counterpart matching rates be- 
tween these rejected sources and optical source catalogs are con- 
sistent with those expected solely from random coincidence, which 
suggests that the majority of these sources are indeed spurious and 
that few valid AGN are removed by this procedure. This evidence is 
strengthened by visual inspection of the rejected sources, which in- 
dependently suggests that the vast majority of the rejected sources 
are indeed spurious. 



4.2.2 Improving Source Positions 

After removal of the most likely spurious candidate sources, 
source and background regions were redetermined for all surviv- 
ing sources. These source regions were analyzed carefully to esti- 
mate the best centroid positions. The source positions determined 
by WAVDETECT do not utilize the full information about the lo- 
cal PSF, in particular at large off-axis angles where the shape of 
the PSF becomes complex. This is especially true for sources with 




o 



Figure 2. An example of spurious candidate sources initially identified 
within the central ~ 1' of MACS J 193 1.8-2634. Contours show the respec- 
tive PSF region (~ 90% EEF) for each candidate source. This particular 
cluste r shows evidence for both a recent merger and powerful AGN feed- 
back (ghlert et al. 2 01 it) , which results in the sharp fronts and cavities in 
the ~ 50kpc surrounding the cluster center. The black source has been 
identified as spurious automatically by ae (Pb ^ 10~ 3 ), while the white 
sources had to be removed manually. Although these sources are signifi- 
cant compared to their local background, they are not AGN. Such candi- 
date sources can only be reliably removed by visual inspection. The only 
true AGN source included in the final catalog is the magenta source at 
the center. This is one of only two luminous cD galaxies in this sample 
(L x ~ 8 x 10 43 ergs-' ). 



"mixed" observations (using both ACIS-I and ACIS-S), where no 
PSF information could be input to WAVDETECT in our analysis. 
We used as to improve the source positions, utilizing three different 
algorithms: (1) calculating the centroid from the data (DATA); (2) 
correlating the PSF within the neighboring region (CORK); and (3) 
reconstructing an image of the region neighboring the source using 
the background and source information and a m aximum likelihood 
algorithm (ML). Tests carried out by the a; team l lBroos et al.l20ld) 
provide guidance as to which of these algorithms is best suited for 
a particular source, based on its position and the positions of neigh- 
boring sources. Following their recommendations, we have moved 
all sources to their DATA positions when they were uncrowded and 
located near the aimpoint (within 8'); to their CORR positions when 
the sources were far from the aimpoint; and to their ML positions 
whenever other sources are sufficiently close that the wings of their 
respective PSF's may overlap. The typical displacement from the 
original WAVDETECT positions is < 1". 

4.2.3 Source Photometry and Spectroscopy 

With the source positions set, final source products were extracted 
from each observation. For each new source position, the PSF was 
simulated by ae at the nominal PSF energy (1.497 keV) as well 
as four other energies (0.277, 4.510, 6.400, and 8.600 in order to 
determine the best estimate of the PSF. Background regions and 
models were recalculated for every source, and spectral products 



including ancillary response files and response matrices were ex- 
tracted, based on the model PSF region. The spectral products, PSF 
models, and no-source probabilities were combined for sources ob- 
served multiple times. In order to determine source fluxes, we uti- 
lized the photon flux measurements for each source computed by 
ajf] and converted them to energy fluxes by assuming a power-law 
model with photon index T = 1.4. This value is consistent with 
the expected photon index for the field population of AGN. Spec- 
tral analyses of the individual sources in our catalog are consistent 
with this assumption. Statistical uncertainties on the source fluxes 
are typically at the ~ 30% level, taking into account uncertainties in 
the source count rates, background, and best-fit power-law model 
parameters. For sources with successful spectral fits, the median 
photon index is T = 1.43, while the mean is T = 1.37. 

4.2.4 Final Source Selection 

Our threshold for inclusion in the final point source catalog is P\, < 
10 r 3 in any band. These thresholds were chosen to be conservative. 
As a comparison, the CDFS uses a threshold of P b < 0.004 in 
any band and led to a sample of high purity (i.e. a catalog with 
few spurious detections expected based on simulations and optical 
counterpart matching rates). A total of 4210 sources are included 
in the final catalog. The distribution of P b in the full band for all 
sources is shown in Figure[T] 

Using estimates of the PSF and local background, we are able 
to efficiently remove the majority of spurious sources identified by 
WAVDETECT, including sources associated with diffuse, extended 
ICM emission. However, we also investigated the central regions of 
every galaxy cluster visually. Despite the best attempts of as to sup- 
press sources associated with the diffuse ICM, structures that vary 
significantly on angular scales comparable to the Chandra PSF can 
be mistaken for point sources at high formal significance, and can 
only be reliably removed with visual inspection. The vast majority 
of these sources are located within ~ 100 kpc of the cluster cen- 
troids, and are pre dominantly located in clusters identified as host- 
ing co ol cores (e.g. lDunn & Fabianl2008ljHlavacek-Larrondo et al.l 
l2012h . Cool core clust ers typically have sharpl y peaked surface 
brightness profiles (e.g. lPeterson & Fab ian 2006). The presence of 
a cool core is also closely linked to AGN feedback by the cluster's 
central galaxy (e.g. Fabian et al. 2003, 2006: McNamara & Nulsen 
20071; ISanders & Fabian 2007; Milli on et al] l20ld : IWerner et al.l 



2010), and ICM substructures arising from this feedback such as 



shock fronts or cavities can also result in structures falsely classi- 
fied as bright point sources. Since the cool cores are typically asso- 
ciated with the centers of clusters where the galaxy density is also 
highest, these sources also commonly have optical counterparts by 
random coincidence. Examples of such false candidate sources are 
shown in Figure[2]fn the extreme cool core cluster MACS J1931.8- 



2634 jEhlert et al 



limb- 

In our visual inspections, we were con- 



servative, removing all candidate sources that could not be clearly 
identified as bright AGN visually. A total of ~ 20 sources were 
removed in this fashion. It is possible, therefore, that our results 

5 We have used both the "FLUX1" and "FLUX2" measurements made by 
as. The "FLUX1" measurement is denned as the net counts in each chan- 
nel divided by the ACIS ancillary response in the channel (for that source 
position in the same energy band) and the exposure time, summed over all 
channels. The "FLUX2" measurement is defined as the net counts in an 
energy band divided by the mean ancillary response of the ACIS detector 
and the exposure time. Both of these measurements convert the measured 
counts into units of photons citT 2 s~ 1 . 



slightly underestimate the AGN density in the innermost regions 
of the clusters (r < 50 kpc). Since the cluster emission is bright- 
est and most structured near the cluster center, the properties of 
all point sources identified in these regions will be subject to sys- 
tematic uncertainties due to potential errors in the background es- 
timation. Based on our conservative visual selection, only two X- 
ray luminous AGN are unambiguously identified within 50 kpc of 
the cluster centroids: one in MACS J 1427.2+4407, and the other 
in MACS J 193 1.8-2634. Additionally, these visual inspections en- 
abled us to manually remove sources that were correlated with in- 
strumental artifacts unrelated to physical point sources. 

4.2.5 Initial Optical Counterpart Matching 

Although we reserve detailed counterpart matching for future 
work, we have performed an initial check comparing the frac- 
tion of matches between our X-ray point sources and cata- 
logs of optical sources der i ved from S u baru i maging in hand 
von der Linden et al] |2012| ; iKellv et al] |2012| ; lApplegate et al] 



I2OI2I) . The optical catalogs include all source types, and are com- 
plete to a magnitude limit of ~ 25 in the Subaru R-haxiA. We 
adopted a constant matching radius of 2" between X-ray and opti- 
cal sources, and in instances where more than one candidate optical 
counterpart within 2" of the X-ray source was present, the brightest 
optical source was chosen as the counterpart. |j In total, the number 
of sources with an optical counterpart brighter than magnitude 25 
{Subaru R-band) is 2826 out of 4210, or 67%. 164 sources in the 
X-ray catalogs reside in regions of the optical images unsuitable 
for reliable photometry (i.e near CCD chip gaps or bright stars), so 
no counterparts for those sources could be established. The coun- 
terpart matching rate between X-ray sources and R-band sources is 
sligh tly lower than th e CDFS at the same flux and magnitude limits 
(e.g. IXueetai1l20T l|), as can be expected given the greater depth 
and wavelength coverage of the follow-up data used for that sur- 
vey. For sources with successful optical counterparts, the median 
separation between the X-ray source and its optical counterpart is 
approximately 0.7". 

4.2.6 Source Extension 

We have performed a check for spatial extension of o ur sources 
using the Kolmogorov-Smirnov (K-S) test discussed bv lXue et al] 
d201ll) . For each detected X-ray source we derived cummulative 
EEFs for both the PSF model and source counts out to 90% EEF. 
We then used a K-S test to compute the probability (Pks) that the 
two sets of cummulative EEFs are consistent with one another. A 
total of 236 sources out of 4210 (5.7%) have a measured value of 
Pks < 0.05 (i.e. have a distribution of source counts inconsistent 
with the simulated PSF at 95% confidence), and a total of 127 
sources (3%) have Pks < 0.01 . These results are consistent with 
measurements from the CDFS dXueetalj|201ll) . We conclude that 

6 We find that a 2" matching radius maximizes the counterpart matching 
rate without introducing a significant number of coincidence matches, and 
is sufficiently large to account for any systematic pointing offsets that may 
exist between the Chandra and Subaru data. 

7 As an example, for sources brighter than 3xl0~ 15 ergcm~ 2 s _1 in the soft 
band, the fraction of X-ray sources with optical counterparts with R < 25 
is ~ 91% (31 out of 34), while ~ 86% of the X-ray point sources in this 
sample at the same flux limit have optical counterparts with R < 25. 



at most ~ 2% of the sources may have spatial extension. Visual in- 
spection of the sources with Pks < 0.01, however, suggests that the 
majority of these sources do not have significant extension beyond 
the PSF, especially for sources observed multiple times at different 
off-axis angles. We therefore treat these results as a conservative 
upper limit to the true number of extended sources in our catalogs. 

4.3 The Final Point Source Catalog 

The final point-source catalog has a total of 42 10 X-ray bright point 
sources detected at our selection criteria (P b < 10~ 3 in any band) 
across a total survey area of 4.2 deg 2 . In the full (0.5 - 8.0 keV), 
soft (0.5 - 2.0 keV), and hard (2.0- 8.0 keV) bands, a total of 4099, 
3344, and 2629 sources, respectively, satisfy the selection criteria 
for that specific band. The majority of these sources (3732) are de- 
tected in at least two bands, and 2019 sources were detected in all 
three bands. Only 256, 81, and 30 sources satisfy the selection cri- 
teria exclusively in the full band, soft band, and hard band, respec- 
tively. Treating the no-source binomial probability formally, there 
is a < 5% chance of there being more than seven spurious sources 
in any band. General information about the point source catalogs 
for each cluster can be found in Table [2] As an example, the ini- 
tial and final point source positions for the galaxy cluster MACS 
J1931.8-2634 are shown in Figure[3] 



5 SENSITIVITY MAPS 

The minimum flux at which point sources can be reliably detected 
in a given Chandra observation and energy band depends on the 
exposure time of the observation, the local background, vignetting 
corrections, the local PSF, the assumed properties of the source, 
and the manner in which point sources are selected. For each clus- 
ter field, we calculate the flux limit at every position in the field 
of view for a canonical AGN source, using the effective exposure 
maps (5j3}, background maps ( i]5.1t . and the no-source binomial 
probability (Equation[TJ. Sensitivity maps were created in all three 
science bands. 

5.1 Background Map Creation 

Background maps for each cluster and energy band were created 
directly from the cluster data, masking all point sources in the final 
catalogs. Masked regions were re-filled using the CIAO routine dm- 
filth, which replaces the excised regions with random values drawn 
from a Poisson distribution, the mean of which is sampled from a 
unique background region for each source. For the excised source 
regions, we used circles with radii 1.5 times the 99% enclosed en- 
ergy fraction radii for each source. For refilling, we used annuli 
with inner and outer radii of 1.5 and 2.5 times the ~ 99% enclosed 
energy fraction radii. Approximately ~ 25% of the pixels in each 
image were re-filled by this procedure. The resulting background 
map is the sum of unresolved cosmic sources, diffuse cluster emis- 
sion, Galactic foreground emission, and instrumental background 

8 Since the no -source binomial prob ability is a formally derived null- 
hypothesis test tWeisskopf et al.|[20Q7t) . the true false-positive rate should 
not deviate significantly from these estimates. This is slightly complicated 
by the possibility of extended sources at high significance, but neverthe- 
less we expect that these catalogs are not contaminated by large numbers of 
spurious sources. 



from each observation. In this study, resolving the relative contri- 
butions of these background components is not important, as the 
sensitivity to point sources only depends on the total background 
present in each pixel. 

5.2 Deriving the Sensitivity Maps 

With the total background maps in hand, we determined the flux 
limit for point sources anywhere in a given cluster field. We deter- 
mine the minimum number of counts required for source detection 
by solving Equation [Tj for our catalog threshold no-source bino- 
mial probability, P b < 10~ 3 . To do so, we require estimates of the 
background counts in both a source aperture (B src ) and an external 
background region (B cxt ), both of which will depend strongly on 
the local PSF. For these calculations, the PSF was assumed to be 
a circular aperture with a radius that depends only on the angular 
distance between the source position and the exposure-weighted 
aimpoint, a p . The radius of the aperture was estimated using the 
~ 90% EEF radius of point sources at the same off-axis angle. De- 
terminations of the background counts in the source aperture, 6 src , 
are taken directly from the background maps discussed above, us- 
ing the local PSF aperture. Estimates of the external background 
counts, B exl , are determined from the B cxt values used in the clus- 
ter point source catalog. We set the B cxt to the maximum B cxt value 
of the cluster-catalog sources that are located in the same annular 
region as the position in question, with the inner/outer radii being 
a p - 1.0' and a p + 1.0'. By choosing the maximum value of B exl , the 
ratio BACKSCAL is maximized for every source aperture, mean- 
ing that the counts required for source detection in an aperture is 
minimized. 

With estimates of fi src and 6 exl for each position in the back- 
ground map, we numerically solve the binomial probability equa- 
tion to determine the minimum number of counts required for de- 
tection under the source-detection criterion. These counts are then 
converted into a count rate given the effective exposure map (g3), 
and the count rate is converted into a physical flux assuming an 
absorbed power-law spectrum with a photon index of T = 1.4. 
The above procedure takes into account PSF broadening with off- 
axis angle, the variations in the effective exposure time (due to vi- 
gnetting and CCD chip gaps, for example), and variations in the 
background across the field of view, including the diffuse galaxy 
cluster emission. An example sensitivity map, for the cluster 
MACS J1931. 8-2634, is shown in Figure g] The impact of the 
diffuse cluster emission on the local flux limit is apparent at the 
center of the sensitivity map, where the local flux limit is a fac- 
tor of ~ 5-10 times higher at the cluster center compared to 
the surrounding region, despite the significantly sharper PSF at the 
cluster center. We define the band-specific flux limit for each clus- 
ter field as the minimum flux in that energy band to which 50% 
of that particular field is sensitive. These flux limits range from 

9 Only one sensitivity map is determined for each cluster, as this procedure 
is equally valid for clusters observed multiple times. The only difference 
is that the background maps and exposure maps from each individual ex- 
posure are summed together before the local flux limit is calculated. The 
median distance between any observation aimpoint and the mean aimpoint 
is ~ 0.7', therefore changes in the local PSF between observations should 
not be significant. For a few clusters, observations have aimpoints appre- 
ciably far away from the mean aimpoint (> 2'), which adds a source of 
systematic uncertainty to the sensitivity maps. These systematic uncertain- 
ties are not large compared to the total survey area over all fields, and are 
therefore inconsequential for the final results. 
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Figure 3. The Chandra full-band counts image for the galaxy cluster MACS J1931. 8-2634 lEhlert et al.l201 ll) . with the initial and final point source catalogs 
overlaid. The magenta ellipses correspond to the initial point sources detected by WAVDETECT, while the blue circles correspond to point sources that satisfy 
the threshold Ft, < 10~ 3 in the full band. The radius of the blue circle corresponds to 1.5 times the 99% EEF radius for the point source position. For this 
cluster, 162 initial candidate sources were detected by WAVDETECT. This was refined to a final sample of 137 high-confidence (Pb < 10 -3 ) sources detected 
across all three bands. 



1-9x10 15 erg cm 2 s 1 in the full band, and are shown for each 
cluster and energy band in Table [2] 



6 RESULTS 

Our final source catalogs and sensitivity maps allow us to deter- 
mine the distribution of X-ray bright point sources across the clus- 



ter fields. We present two main results in this study: the cumulative 
number counts of point sources per unit sky area across the fields; 
and their radial distributions. In both cases, the dominant uncer- 
tainty is the Poisson uncertainties in the total number of sources. 
The expected Poisson fluctuations for a sample of size n is es- 
timated using the l-<x asymmetric confidence limits of iGehreli 




— 




2.40e-16 l^Oe-15 5.02e-15 2.03e-14 

Figure 4. The sensitivity map for MACS J1931. 8-2634 in the 0.5-8.0 keV energy band, calculated by solving for the number of counts required to detect 
a source at a binomial probability of Pb= 10~ 3 . These count estimates were converted to a physical flux assuming an absorbed power-law spectrum with a 
photon index of T = 1.4. For this particular cluster, the minimum full-band flux to which 50% of the area is sensitive is 1.8 X 10~ 15 erg cirT 2 s _I . The dark 
region at the center of the cluster shows that the local flux limit near the cluster is approximately 5-10 times higher than the surrounding region despite being 
the region with the sharpest PSF, due to the presence of the bright cluster "background" emission. 



dl986h . The 1-cr upper confidence limit Au and 1-er lower confi- 
dence limit A L for a sample of n sources are approximated as 



A v (n) = n + 1 + V0.75 + n 



which are both accurate within a few percent for all values of n. 
These Poisson uncertainties, which are at the ~ 5% level, dominate 
over uncertainties in the source flux measurements and uncertain- 
ties in the sensitivity maps. Our results do not strongly depend 

10 This has been confirmed by Monte Carlo simulations, calculating the 
logN-logS number counts from random realizations of the source fluxes and 



Table 2. Summary of AGN source numbers and flux limits for each cluster 
field. The columns are ( 1) cluster name; (2) the number of sources in the ini- 
tial candidate point source catalog for each cluster, produced by WAVDE- 
TECT; (3) the number of sources included in the final catalog followed by 
the number that satisfy Pt, < 10~ 3 in the full band, the soft band, and the 
hard band, respectively; (4) the flux limit for each cluster observation in 
the full, soft, and hard bands, defined as the minimum flux to which 50% of 
the survey area is sensitive, in units of 10~ 15 ergcm~ 2 s -1 . Clusters denoted 
with a dagger (t) utilized a mixture of ACIS-I and ACIS-S observations. 



Cluster Name 


jVwav 


A'cat/Nfull/AWl/A'hard 


Flux Limit 


Abell 209 


155 


56/52/37/34 


5.01/2.00/6.31 


Abell 963+ 


122 


79/77/63/47 


8.91/2.51/14.13 


Abell 2261 


117 


75/72/62/50 


3.98/1.58/5.62 


Abell 2219 


59 


40/35/32/18 


3.16/1.26/4.47 


Abell 2390 


146 


93/93/74/56 


5.62/1.58/10.00 


RXJ2 129.6+0005 


282 


109/104/85/59 


5.01/1.78/7.08 


Abell 521 + 


185 


141/137/113/84 


3.16/1.00/4.47 


Abell 1835 


220 


177/174/139/116 


2.51/0.79/3.16 


Abell 68 


81 


50/49/38/31 


7.08/2.82/10.00 


Abell 697 


156 


69/67/52/43 


5.01/1.78/6.31 


Abell 2537+ 


139 


104/102/81/64 


3.55/1.12/5.01 


MS2137.3-2353 


105 


68/66/54/44 


1.78/0.56/2.82 


MACSJ1931.8-2634 


162 


137/135/105/93 


1.78/0.63/2.51 


MACSJ1 115.8+0129 


109 


94/93/77/63 


2.51/0.89/3.55 


11X11532.9+3021+ 


244 


46/43/38/26 


6.31/2.00/7.08 


ZWCL1953 


145 


106/103/85/52 


3.98/1.41/5.62 


MACSJ0949.8+1708 


96 


49/48/39/34 


5.62/2.24/7.08 


MACSJ1720.2+3536 


150 


125/121/102/83 


2.82/0.89/3.55 


MACSJ1731.6+2252 


136 


83/80/65/49 


4.47/1.58/5.62 


MACSJ22 11.7-0349 


144 


80/77/55/46 


5.01/1.78/7.08 


MACSJ0429.6-0253 


147 


91/89/70/59 


3.98/1.41/5.62 


MACSJ2228.5+2036 


128 


78/76/61/51 


4.47/1.58/6.31 


MACSJ045 1.9+0006 


84 


52/51/40/31 


7.08/2.82/10.00 


MACSJ1206.2-0847 


176 


86/84/59/53 


3.98/1.41/5.01 


MACSJ0417.5-1154 


163 


131/127/102/87 


4.47/1.00/4.47 


MACSJ2243.3-0935 


126 


62/59/49/38 


4.47/1.78/5.62 


MACSJ0329.6-0211 


132 


102/96/77/63 


2.51/0.89/3.55 


RXJ1347.5-1145 t 


366 


86/81/60/57 


2.82/1.00/5.01 


MACSJ1621.3+3810 


186 


159/154/134/112 


2.00/0.71/2.82 


MACSJ1 108.8+0906 


173 


80/79/60/47 


3.55/1.26/5.01 


MACSJ1427.2+4407 


186 


95/95/74/62 


2.82/1.00/3.98 


MACSJ2214.9-1359 


154 


104/104/83/60 


3.16/1.12/4.47 


MACSJ0911.2+1746 


103 


88/85/75/51 


2.82/1.00/3.98 


MACSJ0257. 1-2325 


164 


106/102/85/58 


3.55/1.26/5.01 


MACSJ0454. 1-0300 


169 


85/80/70/48 


6.31/2.24/7.94 


MACSJ1423.8+2404 t 


185 


127/123/105/84 


4.47/1.41/7.08 


MACSJ1 149.5+2223+ 


129 


107/105/85/66 


6.31/2.24/7.94 


MACSJ0717.5+3745 


200 


156/146/133/95 


2.82/1.00/3.98 


MS0015.9+1609 


140 


116/114/94/81 


1.78/0.63/2.51 


MACSJ0025.4-1222 


203 


182/176/150/113 


2.82/0.89/4.47 


MACSJ2129.4-0741 


129 


106/104/88/67 


3.16/1.12/4.47 


MACSJ0647.7+7015 


106 


91/89/68/50 


3.16/1.12/4.47 


MACSJ0744.8+3927 


159 


139/137/111/98 


2.24/0.71/3.16 



on the choice of photon index for the canonical AGN source. In 
all instances, we compare our results to the results from both the 



sensitivity maps. The uncertainties on the source fluxes were taken from our 
spectral fits, while the uncertainties in the sensitivity maps were assumed 
to be ~ 20% and applied coherently across an entire field of view. Both 
effects were shown to contribute negligibly to the overall uncertainty when 
compared to the Poisson fluctuations on the detected sources. 



CDFS and the Cha ndra COSMOS survey (e.g. lElvis et alj[20Qg| : 
iPuccetti et al .120091) in the same energy band. These surveys should 
provide estimates of the number density of field sources not asso- 
ciated with the clusters, at l east within the unce rtainties associated 
with cosmic variance (e.g. iLehmer et al. IHoH. We note that the 
COSMOS results presented here are obtained from re-analyzing the 
Chandra COSMOS field data using our pipeline. This ensures that 
discrepancies between the clusters and field surveys due to differ- 
ences in the analysis procedure and calibration are minimal. We 
omit comparisons to published results based on the XMM-Newton 
COSMOS survey (e.g. lCappelluti et alj2009h due to systematic dif- 
ferences in flux ca libration, which are estimated at the ~ 10% leve l 
in the soft band dNevalainen et al .1 12010| : iTsuiimoto et al.l 1201 ll ). 
Such systematic uncertainties in the source fluxes lead to system- 
atic uncertainties in the cummulative number counts at a given flux 
of approximately ~ 30%. Additionally, Chandra and XMM-Newton 
calibration differences are energy dependent, and will affect the 
hard and soft bands differently. 

In order to ensure survey completeness, we restrict ourselves 
to the survey area within a 15 arcminute radius of the aimpoints 
in calculating the logN-logS cumulative number counts and ra- 
dial profiles. As we will show, we can reliably measure the pro- 
jected source density within ~ 3 r 50 o or 3.5 Mpc of the clusters 
after imposing this cut. The total survey area sensitive to differ- 
ent flux levels after imposing this 15 arcminute cut is shown in 
Figure [5] This shows a steep drop in survey area at fluxes below 
~ 3 x 10~ 15 ergcnr 2 s (0.5 - 8.0 keV). There is comparable sen- 
sitivity to a given flux in the hard and full bands, due to relaxed 
selection criteria in the hard band. 

Only a small number of sources in the final catalog within 15 
arcminutes of the aimpoints have measured full-band fluxes below 
the flux limit at their respective positions (282/4084 ~ 6.9%), and, 
in the majority of these cases the flux measurements are consistent 
with the flux limits within statistical uncertainties. Only 73 (1.8%) 
of the point sources in the final catalog have a measured flux in- 
consistent with the local flux limit at their respective positions at 
a level greater than 68% confidence at their respective positions. It 
is expected that some sources will have such characteristics, given 
differences between the spectra of the point sources and the canon- 
ical AGN spectrum (r = 1.4). The sources with fluxes below their 
local sensitivity limits have spectral fits that are significantly softer 
than the canonical T = 1 .4, with a median value of T ~ 2. 



6.1 Cumulative Number Counts 

The cumulative number density of sources above a given flux (5 ) 
is calculated as 



.(>.)= £ i 



(2) 



where £2, is the total survey area sensitive to the i source flux 
Si. The log TV - log S cumulative number counts for sources in 
the soft, hard, and full energy bands are shown in Figure [6] to- 
gether with the CDFS and COSMOS results in the same energy 
bands dLehmer et al. I l2012h . All fluxes have been corrected for 
Galactic absorption. The cumulative number counts in all three 
bands for the cluster fields, the CDFS, and COSMOS show simi- 
lar shapes, con sistent with the frequently observed broken power- 
law shape (e.g.lCowie et alj|2002l ; (iVIoretti et ai]|2003l ; iBauer et al.l 
l2004lLehmer et al.l2012l)~ A small excess of sources is observed at 
higher fluxes in all three energy bands, and at low fluxes all three 
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Figure 5. Survey solid angle as a function of flux limit in the soft, hard, and 
full bands for all cluster observations. In order to ensure survey complete- 
ness, we have only included the region within 15 arcminutes of the aim- 
points in calculating the survey area presented here. The total survey area 
within 15 arcminutes of the aimpoints over these 43 clusters is 4.06 deg 2 . 



bands are consistent with the expected field population, although 
the largest relative excess is observed in the soft band. 



Our unique knowledge of the overdensity radii, rsoo, allows 
for new insights into the influence of the cluster environment on the 
AGN population. Using Monte Carlo methods, we measure 1.1 ± 
0.6 excess sources per cluster field with respect to COSMOS within 
r soo in the full band, 1.1 ± 0.4 excess sources per cluster detected 
within rsoo in the soft band, ~ 0.5 ± 0.4 in the hard band. 

We have fitted the observed X-ray point source density profiles 
in all three bands with a Power-Law+Constant model. 



N x (r) = No 



(- 



(3) 



We used a Markov Chain Monte Carlo (MCMC) method to deter- 
mine confidence intervals for the power-law index and accounting 
for the source count fluctuations in each bin and the covariance 
between the power-law index and the field density. We assumed a 
uniform prior for the power-law index /? between -5 < yS < 5, and 
sampled over a range of expected field source densities, Cx, using 
a Gaussian prior based on the COSMOS value with a = 10%. This 
value for a is sufficiently large to account for both the statistical 
fluctuations in the COSMOS source counts (~ 5%) and the sys- 
tematic uncertainty due to cosmic variance. The resulting posterior 
distributions in each energy band are statistically consistent with 
one another. For the full, soft and hard bands the resulting posterior 
distributions measure the power-law index as p = -0.70 ± 0.23, 
P = -0.52 ±0.34, and/8 = -0.75 ±0.41, respectively. The results of 
fits to all three bands, the expectations from the field surveys, and 
determinations of the number of excess sources are shown in Table 

m 



6.2 The Radial Distribution of X-ray Sources 

Measurements of the projected source density as a function of 
cluster-centric radius provide both a powerful investigative tool for 
understanding the causes of AGN activity and an important system- 
atic cross-check to ensure that detections of excess AGN sources in 
the cluster fields are not due to fluctuations in the field populations. 
The radial distribution of point sources about the cluster centers has 
been calculated for all point sources with full-band fluxes above 
5 x 10~ 15 ergcirr 2 s" 1 . Similar analyses were performed in the soft 



band and hard bands, with flux limits of 3 x 10 15 erg cm 2 s 1 



and 



10~ 14 ergcirr 2 s" 1 , respectively. For all of these flux limits, there 
is no evidence for survey incompleteness out to 3r 50 o- A total of 
2056, 968, and 930 sources satisfy these criteria in the full, soft, 
and hard bands, respectively. The radial distributions are plotted in 
two ways: firstly, in terms of the metric distance (r, in Mpc, Figure 
|7(a)fr from the cluster center; and secondly (and more usefully), in 
terms of the overdensity radius (r, in rsoo, Figure [7(b)) . The radii (in 
Mpc or r 5 oo) of each source were calculated assuming the sources 
are at the same redshift as the cluster. In both representations, there 
is evidence of an excess of point sources in the central regions of 
the clusters. The radial profiles for the soft band and hard band are 
shown in Figures |8(b)| and |9(b)| respectively. 

At large radii, the measured source number densities con- 
verge to a constant source density, consistent with the measured 
value from CDFS and COSMOS. Within the virialized cluster re- 
gion (rsoo) an excess of sources is observed in all three bands. We 
also observe tentative evidence for a second peak in the profiles be- 
tween ~ 1.5-2 rsoo, especially in the soft band profile. The overall 
excess of sources above the field population observed in the hard 
band is smaller than in the soft and full bands, as expected given 
our hard band flux limit. 



6.3 Comparisons with the Cluster Galaxy Distributions 

Since galaxy clusters are, by definition, regions of high galaxy den- 
sity, it is likely that any excess of X-ray point sources observed in 
cluster fields will, at least in part, simply be a result of the overden- 
sity in the underlying galaxy population, rather than influences of 
the ICM on the evolution of galaxies. Formal counterpart matching 
between optical a nd X-ray sources in these clusters, such as thos e 
performed in (e.g. lMartini et al]200ll2009llGalametz et alj2009h , 
will be left to future work, but we can nevertheless compare the 
radial distribution of X-ray sources in our fields with the expected 
radial distribution of cluster galaxies. Although the AGN fraction 
cannot be directly determined with such a comparison, we can de- 
fine an effective AGN fraction, d> as 



Nx(r) - C x 
N (r) 



(4) 



where No(r) is a model describing the projected galaxy density, 
N x (r) is the number density of X-ray point sources and Cx is the 
background (field) number density. This function should be repre- 
sentative of general trends in the AGN fraction profile, albeit with 
arbitrary normalization. We perform our analysis within our frame- 
work discussed in i]6.2| to fully account for the covariance between 
m odel parameters. We adopt the shape for N (r) from the results 
of lPopesso et al.| j2007t). who identified cluster member galaxies for 
217 X-ray selected clusters using observations from the Sloan Digi- 
tal Sky Survey binned b y mass, using a magni tude cut of r < -18.5. 
The clusters included in lPopesso et al]( l2007l) span a large range of 
cluster masses, and like the clusters in this sample were selected 
based on their X-ray properties. Those authors use a King model 
parameterization to model the projected galaxy distribution. For the 
most massive galaxy clusters in their sample (Js 5x 10 I4 M G , com- 
parable to the masses of clusters in this sample), their model has a 
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Figure 6. Cumulative number counts (logN-logS) in all t hree energy bands fo r the cluster field X-ray point source catalog in black. In red are the cumulative 
number counts fo r the CDFS in the same energy band iLehmer et al. 2 0121) . and in blue are the results of our re-analysis of the Chandra COSMOS field 
lElvis et al. 2009). We only include bright sources in determining the radial distribution of AGN, the band-specific flux limit of which is denoted by the 
vertical dashed line in each figure, (a): The full-band cumulative number counts, (b): The soft-band cumulative number counts, (c): The hard-band cumulative 
number counts. All three of these logN-logS curves show similar trends, with a small excess of sources in the cluster fields detected at high fluxes. 



Table 3. Choices of flux limits and expectations from control fields for the radial profiles in all three energy bands. The columns are: (1) The energy band; (2) 
the flux limit imposed in the radial profile fits, in units of ergcm~ 2 s -1 ; (3) the density of field sources from the CDFS at that flux limit; (4) the density of field 
sources from COSMOS at that flux limit; (5) the number of sources detected within /-500 at that flux limit, across all clusters; (6) the survey area within r^QQ at 
that flux limit, in units of deg~~; (7) the average excess number of sources per cluster, determined by extrapolating our measurements of the field density from 
COSMOS to each radial bin within rsoo; and (8) the best-fit power-law index for the projected source density of cluster member AGN. 



(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


Band 


Flux Limit ( erg 


; cm 2 s 1 ) 


CDFS(deg~ 2 ) 


COSMOS (deg -2 ) 


"500 


n 50 o (deg -2 ) 


Excess 


P 


Full 


5x irr 


15 


690 ± 75 


672 ± 29 


565 


0.77 


1.1 ±0.6 


-0.70 ± 0.23 


Soft 


3x irr 


15 


250 ± 45 


255 ± 18 


244 


0.78 


1.1 ±0.4 


-0.52 ±0.34 


Hard 


1 x irr 


14 


220 ± 45 


287 ± 19 


245 


0.78 


0.5 ±0.4 


-0.75 ± 0.41 



core radius of r c = 0.24 ± 0.02 rw 1 .['"lThey also fit their proj ected 
galaxy data with a projected NFW dNavarro et al J 19951 , 1 1997b pro- 
file with a concentration parameter c = r 2 m/r s = 4.2 ± 0.3. This 
result is similar to o ther studies that use X-ray selected samples of 
galaxy clusters (e.g. lCarlberg et al.l 19971) . We marginalize over the 
full range of acceptable King models. Within the framework of this 
analysis, we divide each realization of the power-law models de- 
scribing the X-ray sources by these two models for the projected 
galaxy density, and determine the full posterior distribution for (f> 
for each radial bin. The median value of the posterior distribution 
at each radial bin in the full band, along with its corresponding 68% 
confidence interval is plotted in Figure[l0] The resulting profile for 
(f>(r) is shown to be clearly increasing monotonically with radius. 
Compared to the value of <f> at ~ 0.25r 50 o, the data suggest that the 
AGN fraction is roughly 2-4 times larger at 7-500- A similar trend 
is observed when this same test is performed with the projected 
NFW models described above. The data presented here do not 
require any excess sources out beyond r 50 o, but a larger sample 
of galaxy clusters with more additional statistical constraining 
power will be able to determine whether or not this trend con- 
tinues out to and beyond the cluster virial radius. 



7 DISCUSSION 

We have demonstrated that an excess of X-ray point sources is ob- 
served in the directions of 43 massive clusters with Chandra, when 
compared to expectations from deep and medium-deep fields. Al- 
though excesses of X-ray point sources in cluster fields have been 
reported in the literature, our analysis has revealed some details not 
previously discussed. 

The mean excess of sources in these clusters within rsoo is 
~ 1 per cluster field . This is consistent with the results of some 
previous studies (e.g. iGilmour et al.l2009h . These point sources are 
very likely AGN associated with the clusters, with a typical X-ray 
luminosity in the full band (L x (0.5 - 8.0 keV)) of ~ 10 42 ~ 43 erg s" 1 . 

The projected X-ray AGN source density is highest at the cen- 
ters of clusters, with excess sources observed out to ~rsoo- The sim- 
ilar excess of sources seen in the soft and full bands given our 
imposed flux limits (3 x 10~ 15 ergcnr 2 s -1 in the soft band and 
5 x 10~ 15 erg cm 2 s -1 in the full band) suggests that the cluster 



11 We assume that r2oo/ ^500= 1-5, which is a good estimate for a large 
range of co ncentration param eters for dark matter halos with an NFW den- 
sity profile dYang et al.l2009l) . 



AGN may have comparable levels of intrinsic obscuration as field 
AGN. 

The excess sources associated with the cluster are more 
broadly distribute d than the galaxy di stributions in typical X-ray 
luminous clusters jPopesso et al .II2007T) . Our initial modeling sug- 
gests that, within ~ r 50 o, the distribution of AGN can be well de- 
scribed by a power-law with an index of j} 0.7. More detailed 

analysis will require a larger sample and/or improved rejection 
of non-cluster mem ber X-ray sources. Ho wever, we caution that, 
while the clusters in IPopesso et alj d2007l) do span a similar mass 
range to the presen t sample, they are at lower redshift (mean red- 
shift of z ~ 0.1 for lPopessoetall d2007t) versus 0.2 < z < 0.7 for 
the clusters presented here). Blue galaxies are also generally more 
broadly distributed than red galaxies. Nonetheless, our initial re- 
sults suggest that the intrinsic X-ray AGN population in clusters 
has a broader spatial distribution th a n both the red and blue cluster 
member galaxies in IPopesso et al. I J2007h f T2 l A power-law model 
of P ~ -0.7 is broader than nearly all King or projected NFW 
models with a core/scale radius less than ~ rsoo- Future work will 
explore the properties of the optical counterparts for these X-ray 
AGN using the Subaru imaging in hand, providing insights into the 
host-galaxy properties, and allow us to determine the AGN fraction 
more robustly. 

We see tentative evidence for a second peak (~ 2a signifi- 
cance above the COSMOS density in the soft band) of sources 
near the virial radii of the clusters, between 1 and 2 r 5m . Further 
analysis on a larger sample of clusters will be required to con- 
firm the presence of this second peak, given the expected statistical 
fluctuations for each radial bin. A similar feature was reported by 
iRuderman & Ebelind d2005l) . who utilized soft band (0.5-2.0 keV) 
measurements and a cluster sample that overlaps with those in- 
cluded here. This fe ature was not observ ed in the large cluster 
sample analysed by IGilmour et al. I d2009h , who utilized the full 
(0.5- 8.0 keV) band. 

Taking these results at face value, the rising AGN fraction with 
radius suggests that the evolution of X-ray AGN in galaxy clusters 
is qualitatively similar to the evolution of of star-forming galax- 
ies. An evolution of the average X-ray AGN fraction in clusters 
with redshift has been cl aimed based on comparisons between high 
and lo w redshift clusters dMartini et al.l2007tT2009l ; lGalametz et alj 
l2009h . Our study reveals an effect correlated with clus ter radius, 
which in turn relates to the time since galaxy infall dGao et all 
120041) . Importantly, this trend is qualitatively similar to trends in the 

12 X-ray AGN in t he field are most commonly observed in more massive 
redder galaxies (e.g. Xue et alj201fj|) . 
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Figure 7. The projected density of X-ray point sources in the full band (Fx(0.5 - 8.0 keV) > 5 X 10~ 15 erg cm -2 s _1 ), in units of deg~ 2 . For both figures, the 
solid red line denotes the best fit source density of COSMOS at the corresponding flux limit, while the dashed red lines correspond to its 1-cr confidence 
limit, (a): The projected density of sources as a function of radius in units of Mpc, assuming each point source is at the same redshift as the cluster. A total of 
2056 sources were included in the calculation of this profile. A large excess is seen in the central ~ 250 kpc of the cluster centers, (b): The same surface density 
as (a), instead normalized to each cluster's characteristic radius r^QQ. In these units the central density peak is more broadly distributed throughout the central 
~ 0.5 r5oo of the clusters. In both cases, the measured density of background sources is consistent with expectations from COSMOS. The black curve shows 
the most probable projected source density profile from our MCMC analysis; a power law with index fS = -0.7. The background source density for 
this model is slightly lower than but statistically consistent with COSMOS. 
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Figure 8. The projected density of X-ray bright point sources in the soft band (Fx(0.5 - 2.0 keV) > 3 X 10~ 15 erg cirT 2 s ). For both figures, the solid red 
line denotes the best fit source density of COSMOS at the corresponding flux limit, while the dashed red lines correspond to its 1 - <x confidence limit, (a): 
The surface density of X-ray bright soft band sources as a function of radius, in units of Mpc. A total of 968 sources were included in the calculation of this 
profile, (b): The surface density of X-ray bright soft band sources as a function of radius, in units of r^QQ. Just as in the full band, using the overdensity radius 
f - 5oo results in a broader central excess of sources. In both profiles shown here, there is evidence for a second excess near the cluster virial radius in addition 
to the central peak. The background source density is consistent with expectations from COSMOS. As in the full band case, the most probable model from 
our MCMC analysis is overplotted as the black curve (J3 = -0.52), and the background source density is statistically consistent with COSMOS. 
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Figure 9. The projected density of X-ray bright point sources in the soft band (Fx(2.0 - 8.0 keV) > 10 erg cm s ). For both figures, the solid red line 
denotes the best fit source density of COSMOS at the corresponding flux limit, while the dashed red lines correspond to its 1 - <x confidence limit, (a): The 
surface density of X-ray bright hard band sources as a function of radius, in units of Mpc. A total of 930 sources were included in the calculation of this profile, 
(b): The surface density of X-ray bright hard band sources as a function of radius, in units of r^gg. Just as in the full and soft bands, using the overdensity 
radius rsoo results in a broader central excess of sources. The background source density is consistent with expectations from COSMOS in both profiles. 
As in the full band case, the most probable model from our MCMC analysis is overplotted as the black curve (ft = -0.75), and the background source 
density is statistically consistent with COSMOS. 



fract ion of both optically luminous AGN and star-forming galaxies 
(e.g. lLarson et alii 1980fc iBekki et al.ll2002l; IWeinmann et al JlioOrl 



120091 : Ivan den Bosch et alj l2008jlbl : Ivon der Linden et al.l I2010T) . 



Although the data presented here do not measure the absolute frac- 
tion of cluster galaxies hosting X-ray AGN, they suggest that the 
ICM may suppress X-ray AGN activity as galaxies are gradually 
stripped of their fuel reservoirs during infall. Such a result is con- 
sistent with a scenario where the X-ray emission of AGN in clusters 
is linked to the accretion of cold gas. 

Our results indicate that the central regions of clusters are rel- 
atively inefficient at producing X-ray AGN, in contrast to the pres- 
ence of radio bright AGN at the centers of many of the clust ers in 
this s a mple as well as g i ant elliptical galaxies in general (e.g . Best! 



20041; iBest etafl 



Best & Heckmanl 



II as giant elliptical galaxies in general (e.g. aesl 
2007; iDunn & Fabian] 120081 ; iDunn et alj|20ld : 



20121) . Other studies have shown that active cD 



galax i es are c ommonly ra d iatively inefficient ( e.g.lDi M atteo et al 



200 ll |2003|; iBestl |2004 lAllen et al.l |2006| ; iTavlor et alj 1200 
Best et al. 1 120071) . with radiative emission accounting for ~ 1% or 



less of the total AGN energy budget. 

Further analyses utilizing a larger sample of X-ray clusters 
and detailed optical follow-up are underway. Those will provide 
improved measurements of the radial source density profile, and 
sufficient data to determine whether the X-ray AGN population in 
clusters evolves with redshift. Counterpart matching between X-ray 
sources and optically selected galaxies will be needed to measure 
the absolute fraction of galaxies hosting X-ray AGN directly and 
identify which X-ray AGN are cluster members. Measurements of 
X-ray spectral properties such as the average photon index or X- 
ray/optical flux rati o offer important clues as to the accretion pro - 
cess for AGN (e.g. lShemmer et al.ll2005l 120081 : IScott et alj|201ll) . 
and probe the impact of the ICM on the accretion process of AGN. 
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Figure 10. The effective AGN fraction as a function of radius, <Kr),in the 
full band (F x (0.5 - 8.0keV) > 5 X 10~ 15 erg cm" 2 s" 1 ). The dashed line 
corresponds to the median value of <j> in that radial bin, with the solid lines 
corresponding to the 68% confidence interval. We calculate <f>(r) using an 
MCMC analysis that determines the most likely power-law model param- 
eters for projected X-ray source density, taking into account all uncertain- 
ties in and covariance with the density of field sources. These power-law 
models were th en divided by an ass umed model for the projected surface 
galaxy density (Popesso et al. 2007); a King model with a core radius of 
r c = 0.24 ± 0.02 r^og. The general trends suggest that the AGN fraction 
at rsoo is approximately ~ 2 - 4 times higher than the AGN fraction at 
~ 0.25r 500 . 



8 CONCLUSIONS 

We have constructed rigorous catalogs of X-ray bright point 
sources within the fields of 43 massive, X-ray luminous galaxy 
clusters. We have compared the observed X-ray source population 
to deep field surveys. Our main results can be summarized as fol- 
lows: 

1) We measure an excess of point sources detected in the clus- 
ter fields when compared with expectations from deep field sur- 
veys in the soft (0.5 - 2.0 keV), hard (2.0 - 8.0 keV), and full bands 
(0.5 - 8.0 keV). On average, this excess amounts to approximately 
1 source per cluster field. 

2) The excess AGN are primarily located within ~ rsoo of the 
cluster centers. At large radii, the point source density is in good 
agreement with expectations from the CDFS and COSMOS in all 
three bands. 

3) The radial profile for the projected source density can be 
described by a power-law plus constant model, with f3 0.7, con- 
sistently measured across all three energy bands. 

4) A preliminary comparison of our measured AGN profiles 
with the typical distributions of m ember galaxies in co mparably 
luminous X-ray selected clusters l lPopesso et al.l [20071) suggests 
that the fraction of cluster member galaxies hosting X-ray AGN 
rises with increasing radius. At face value, this suggests that the 
ICM may suppress X-ray AGN activity as galaxies are gradually 
stripped of their gas reservoirs during infall, qualitatively similar 
to the strang ulation processes observed f or star-forming galaxies in 
clusters (e.g. lvon der Linden et alj2010l) . 

Further studies are underway, employing optical counterpart 
matching and follow-up spectroscopy. This should enable us to in- 
vestigate the optical properties of the host galaxies, the X-ray spec- 
tral properties of the AGN, and provide new insights into the fuel- 
ing and evolution of X-ray bright AGN in clusters. The full catalog 
of X-ray sources will be made public at a later date. 
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